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deJ™ino jeCtiVe -° f this Study is to experimentally 

an empirical model of the vibrational dynamics of 
far-i l ?f° eCr mK t control Laboratory Experiment (SCOLE) 
facility. The first two flexible modes of this test artirip 
are identified using a linear least-square identification 

ob?a?n U H e H and the . data utilized for this procedure ^re 10 " 
by exciting the structure from a quiescent state 
with torque wheels. The time history data of rate gyro 

sxSua?ion d ^ C ?: ler 0 m f? rS due t0 Citation a£d 
xcitation in terms of free-decay are used in the parameter 
estimation of the vibrational model. P 

The free-decay portion of the data is analyzed usina the 

order^to -^ermine the o^ttaa? Se! 

, . use _ in modelling the response. Linear least-sauarp 

the outDut S 0 f h an A S ? d t0 selact the parameters that best fit 
the output of an Autoregressive (AR) .model to the data ThP 

contro! effectiveness of the torque wheels is thin 
determined using the excitation portion of the test data 
again using linear least squares. ' 
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INTRODUCTION 


Future NASA space mi ?=i 0 £ a “ a ^"?ttctt? struStwal 
highly flexible , SSSSSrls would require 

ver^stringent these 

reqmremen . designed with very accurate 

objectives will nave co ue .V. +-»,*» structure to ensure 

knowledge of the and 

performance robustness to a var Y engineers that 

uncertainties. It is recogni^eu * * x. he d esiqn of active 
there are certain inherent Proems inthj^de J^ ecraf t . 

controllers for this clas <3 desire to offer a means 

Because of these concerns and of the desire to - 

oftomparing technical approaches directly a technics? 

Design Challenge m was ^| c °“| r control Bran ch at NASA 

?he C NAlA/ IEEE Design challenge to C oSteS“522ign 

tesl article? lessors, 

spacetntenna^attached tothetpace Shuttle orbiter by a 

St a°multi- inputtutputttructural dynamics' syste^can 

Tndtot: communication^capability. Much 

?= n : r - a5 rs- I s 1 sjsrs sssa .« 

S^ive workshops held since the conception of SCOLE in 
1983 . 


SCOLE APPARATUS 

The SCOLE hardware and software support is ^scribed in 

affAss-iStSSSsasS.* 

is fixed to ground and only the mast and retiBc^oT pu 
are dealt with in this experiment. The system actuators 

rodite ^r^etntK;r^aty d o^onal“S?? tga. T he 

mounted S rotational ^X^ofand^lh^sl-^ited a'nd 
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cabled 5 no 5 £ 3 ™«- l ! ke operating system version 

of c and Fortran 0 ?? 3 ™ 1 " 9 1S accom Plished in a combination 
diaital Programs. The computer has analog-to- 

data S diotta! ?? ers / used for sampling the rate sensor 
i-™i< dlg i tal to_anal °g (°/ A ) converters used to command the 
orecicp 0 torque w hee Is, and a process timer which achieves 
precise internal timing of the data sampling process. 

SYSTEM MODELS 

The model we seek for SCOLE should incorporate the actual 
natura! frequencies, damping ratios, and control 
e fectiveness coefficients of the system. The viscous 

th^i^modl 1 bS ,? 0d **? d in tSrmS ° f r the daitl P in< ? ratio of 
the 1 mode. To this end, each mode of the vibrational 

dynamics of SCOLE is modeled as a single-input, single- 
output system [3] described by the state-space equation 

£(t) =4x(t) +Bu(t) , (!) 

where 


&= [n^] T 


0 1 
to 2 -2C&) 


P-0 

B ~f 


In this equation, 


co 2 = — 


2Co> = - 


where 


x - modal state vector 

u - control input of reflector end reaction torque wheels 

b control effectiveness parameter of actuator location 
w - natural frequency of mode location 
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f — damping ratio of mode. 

The output is of the form 

y=Hx 


for a rate sensor, 


H= [0, c] , 


( 8 ) 


where c is the mode slope at the sens ° r tai 

To obtain a f Lsomed lonstant over 

the P computlr sample time intervalofTsecondsand^the 

continuous-time model is converted to its discrete tim 

e< ^ U Thus? n the^difference°equation^describing* the motion 

appears as 


x £±1 =4>x r +ru jr . 


( 9 ) 


where 


- a AT 


<j) = e 


( 10 ) 


r=f e At dtB= (4 >-I) 


( 11 ) 


since A is nonsingular^and^is ( the ^^^"tity^atrix. 


ice id 

ThP model used in this work is a ixi^u M r7‘T , ^ 

roefficient difference equation • To accomplish such 
coefficient, .."™‘ essive 4 forl of the discrete time mod 

found bv taking the z-transform of the last equation an 

solving for the sampled sensor output, y^, in terms 

input actuator. The auto-regressive moving average (ARMA) 

model appears as 


where 


<a 2^Jf-2 + ^ b l U K-l +b 2 U K-2^ b TW- 

(12) 

a i = 4 > H +< t ) 22 

(13) 

a 2 =< t > 12^ ) 21 _( t > ll < ^22 

(14) 

i?i=r 12 

(15) 
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•^2 $ 12^22 ^ 22^12 


(16) 


b m =c‘b. 


( 17 ) 


■ Ide ]?u lf Ration of the ARMA model parameters is performed 
using the linear least-square estimation (LSE) algorithm. 

method was selected because of its computation 
efficiency and implementation simplicity. The error 
equation used in the estimation is defined to be 


e K=y K - [ a&K-l + * 2^-2 + ( h X Ujr-1 +b 2 U K _ 2 ) b m ] . 


( 18 ) 


The sum of the squared error, 


J= 


i N 


K-0 


'K‘ 


( 19 ) 


is the performance measure to be minimized with 
the parameter desired to be found. 


respect to 


PARAMETER IDENTIFICATION TECHNIQUES 

.. L j^ r least-square estimation (LSE) is used to identify 
the ARMA model parameters. This method is selected because 

The efflclen ?y and implementation simplicity. 

The identification process is carried out for each reaction 
torque wheel and for each mode. The test data in the 

i ?^ 1 fi Catl ° n p f? cess of the SCOLE problem is processed in 
a two-step operation. 

,p M j 0r th e first step, the AR coefficients a, and a 2 of the 
a ^ e lde ntified using the free-decay portion of the 

ooit?nn e nf ?£% ?he s pectral content of the free-decay 
portion of the data is examined using the Discrete Fourier 
Transform and the Hamming window. The free-decay data are 

no ^ se ? nd si 9 na ls due to any modes not 
t 5 e 1 mod f 1 * The filtered data is processed using 
the stantiard least-square estimation to identify the AR 
efficients, a, and a 2 , of the modelled mode. The 
identification of the a, and a 2 coefficients for each mode 
generally depends on the data base used in the estimation. 

conv ^ 6 s ^ gniflcant data are added, the estimates should 
converge to a value and the variance of the estimates will 

ancTthe model Whlch depends on the measurement noise 

and the model. After convergence occurs, the mean and 
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iiSIJgilMgfe" 

control effectiveness coefficient b™ of the torque wn 
n . p(1 is determined. Once the values of the a t and a^ 

estimate of b w is taken into consideration. 


EXPERIMENTAL PROCEDURES 

Experimental work and testing can be conducted on SCOLE 

%; e 

S«sint^ a ii ?«S™o= ?fit nKsA and 

tr ° m A manual C S tructural SgKiS teltls pelformed wherein 
data is collected and analyzed to verify physical modal 

directions and £^£^£-£"3 ^ initial sensors 
simulations of SCOLE. Bias reauiny .. h reference 

are always taken before each run to establish a rereren 

h&httXSlflk £ ££?3 - - a a^-Lis 

o^untilHthe^isplacement^ngle^about'the^-axi^reaches 

?Ive degrees. When the reflector is released, free decay 
da^a recorded and collected. The sane process occurs for 
d a } %oat?no with the reflector being released from the 
mod . f , h ,^' v an d +v-axis directions as shown in Figure 

f" Datfis^onecSd Ld X ierified against predicted natural 

tre ?r$ib?afionIf lyna^ic midel we see* for SCOLE 
incorporating . the actual natural /^uencies damping^^^ 

in^this r k C brusing f th“end mounted reaction wheels 

to excite the structure. Structural excitation tests are 
individually conducted wherein the structure, lnitia y 
rest V issinusoidally forced by a single reaction wheel for 
30 siconds at the predicted mode of interest. The data 
recording is continued for 60 seconds to obtain free decay 
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data. Tests are carried out for each of the x,y, and z-axis 
in ?able ? f ° r b ° th modes 1 and 2 and are summarized 


DISCUSSIONS 

Experimental data and results of the system 
identification process are summarized and illustrated in 

J ab * e l Flgures 5 through 10. Tests and analysis using 
each of the x, y, and z-axis torque wheels at both the first 
and second modes, .4401 and .4764 Hz respectively are 
carried out The identified parameters a 1( a 2/ bj, and the 
computed values of £ and f for all of the tests i^-e 
tabulated in Table 2. 

Figure 5 shows the input excitation signal of the 

whee1 -. Before each test, the structure was 
steadied and bias readings of the sensors were taken and 
accounted for. A sinusoidal forcing signal of amplitude 20 
at the desired mode test frequency was applied to the 
structure for 30 seconds. For mode 1 frequency of .4401 Hz 
the x-axis reaction torque wheel excited the structure the ' 
most. This was in agreement with the predicted first 
bending mode shape shown to occur closer to the pitch or v- 
a ** s dlrect ion. , The y-axis torque wheel had the greatest 
effect on exciting the structure for mode 2 tests, which 
agreed with the predicted second bending mode shape 
occurring in the direction of the roll or x-axis. The z- 
axis torque wheel had the least effect on exciting the 

or ^ bon t- S t-h aS ^ ten< ^ ed to . excite the structure in the yaw 
or about the z-axis direction. Since the first two 

vibrational modes occur dominantly in the x-y plane the z- 
ax 1S torque wheel had little effect on exciting the 
structure at the first two modal frequencies. 

Figures 6 and 7 illustrate the frequency spectral 
ysis on the free decay portion of the data using the DFT 
and Hamming window. The structure would non-periodically 
vibrate in both translational and rotational directions 

about 3 ° seconds of free decay. Therefore, only the 

were considered due to other modes becoming 

bv Jho Li" decay * The length of the data record used 
^ C n t0 be equal to an integer number of 
periods of the sequence. The presence of leakage or 

significant non-zero frequency components occurred when the 
data record was improperly truncated. Therefore, an inteqer 
number of periods represented by the value of N was chosen 
The frequency magnitudes were also proportional to the 
number of periods included in the record length. The 

?r?n^ the n r ber ° f periods ' the larger the magnitude, 
introdnrp^ £ ed ? Ce the effects of the discontinuities 

troduced by truncating the sequence, Hamming windows were 
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used. The Hamming windows represent a noticeable d , bes 
:"f?hrunLnSernSn-z«o^requen?iL\hi?e broadening the 

" in detecting^and^identifying^th^first 

and second modal frequencies, which were very clo 

t0 ^Fiqures 8 and 9 of each test provide a description of the 
parameter i dentif ^^ 

^e^i^orofeL“SrSr!Lr?hr:eireratir£?rom?heir 

arithmetic average with a convergence ^ "* d *\°| t 5 are 

KoSrln TaMeT" S^tEa? frequency and its damping 

S“o c =?n b ?ib?e u r s: S gnft| 3 ' 

identif icatio^techniqueT ^he^SqSency calculated from^he 

show^o^the 1 ^ igures^^he^characteristics^and accuracies of 

the LSE algorithm and the DFT accounted for a part 

dit F?gu?fio of each test shows the control =H™ss 
parameter of the reaction torque 

small when compared “ ^'^fJ^th^diffSence between the 
forced^and'resonant frequencies of "^^^“"co^eTget 

SSffi 

k taken over a converged portion of the data a 

f igures^°*The pSTse^Yf"! dSurbanc^ & TnduceTwhen the 

ulf-r^W^on- e el6Ve 

reaction torque wheel assembly. , T1 ?® r ®^ r ®: t 3 h b a Je for 
seconds of data were not included in the data base 

estimation. 


SUMMARY 

In • thlS f W ?he' fI?sr?io 1 ftex?We 1 modes t ’of V Sc5S i Sas 1 found 
S^ng^L^ineL' Sfst squirl identification procedure. The 

experimentai apparatus and procedures followed and the 
system model assumed were also discussed in this st y. 
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Testing wasd° ne by exciting the structure from a quiescent 
state with torque wheels and recording the time history data 
wfr^Jhp 97 ^ f e " sors and accelerometers. The torque wheels 

Tnn w down and free - d ecay data recorded. The DPT 

and Hamming window were used to analyze the free decay 

portion of the data. The coefficients of an autoregressive 
model to the data were determined using linear least square 
analysis. Next, the control effectiveness of the torque 

da?a aaaL f utn i USi "V he excltati °" portion of the Sst 
da£a'and L llzlng ^ near least squares. Experimental 
data and graphs were also presented to provide a description 

SL dl ? ltal f^ gna; |- Process algorithms and techniques P used 
rmining the first two flexible modes of SCOLE The 
experimental results presented in this paper have th4 
potential to be extremely useful in modelling vibrational 
dynamics of large flexible spacecraft structures. 
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The SCOLE Experimental Apparatus 
Figure 2 


Original photograph unavailable at time of publication 
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Actuator Input 

Figure 5 
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Discrete Fourier Transform 
Figure 6 
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Hamming Window 
Figure 7 
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VARIANCE A1 • .0001020 
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Figure 8 
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Arma Coefficient A 2 

Figure 9 
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MEAN B • 1.5270996 
VARIANCE 8 • .1837367 
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TIME (SECS) 

Control Effectiveness Parameter B 
Figure 10 
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Table 1 


Test Parameters 


TEST 

4 

REACTION 

TORQUE 

WHEEL 

MODE 

EXCITATION 

TIME 

(SECS) 

FREQUENCY 

(Hz) 

amplitude 

FREE 

DECAY TIME 
(SECS) 

SAMPLE 

INTERVAL 

(SECS) 

1 

X-AXIS 

I 

30 

0.4401 

20 

60 

0 1 

2 

Y-AXIS 

1 

30 

0.4401 

20 

60 

0.1 

3 

Z-AX1S ] 

l 

30 

04401 

20 

60 

0.1 

4 

X-AXIS 

2 

30 

0.4764 

20 

60 

0.1 

5 

Y-AXIS 

2 

30 

0.4764 

20 

60 

0.1 

6 

Z-AXIS 

2 

30 

0.4764 

20 

60 

0.1 
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Table 2 


Experimental Identification Results 


■ 

n 

MODE 

■ 

a 2 

GRAPH 

FREQ 

(Hz) 

LSE 

FREQ 

(Hz) 

c 

b TW 

l 

X-AXIS 

] 

I 5235475 

-0 6120654 

0 44715 

0.5351800 

0 7299576 

-5.5898004 

2 

Y-AXIS 

1 

1.6931082 

-0.7789838 

0 44715 

0.4978276 

0.3992480 

0.0080636 



1 

1 5868783 

-0.675171 1 

044715 

0.5223736 

0.5983680 

-0.8804742 

■ 

X-AXIS 

2 

1 8732456 

-0.6641920 

0.46218 

0.5164785 

0.4766312 

0.6327561 

5 

Y-AXIS 

2 

1.7073690 

-0.7947280 

0.46218 

0.4997192 

0.3658724 

1.5270996 

6 

Z-AXIS 

2 

1.6953652 

-0.7836996 

0 46610 

0.5042274 

0.3846554 

-1.4330764 


259 



























































